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ABSTRACT 

We report the detection of four new extra-solar planets from the Anglo-Australian Planet Search 
orbiting the somewhat metal-enriched stars HD 73526, HD 76700, HD 30177 and HD 2039. The planetary 
companion of HD 76700 has a circular orbit -with a period of 3.98 d. With Msini=0.197±0.017Mjup, or 
0.69 times the mass of Saturn, is one of the lo"west minimum mass extra-solar planets yet detected. The 
remaining planets all have elliptical orbits "with periods ranging from 190.5 d to 4.4 yr. All four planets 
have been found orbiting stars from a sub-sample of t-wenty metal-enriched and faint (V<9) stars, "which 
"was added to the Anglo- Australian Planet Search's magnitude-limited V<7.5 main sample in October 
1998. These stars "were selected to be metal-enriched on the basis of their Stromgren photometry, and 
their enrichment has been subsequently confirmed by detailed spectroscopic analysis. 

Subject headings: planetary systems - stars: individual (HD2039, HD73526, HD30177, HD76700) 



small sub-sample of t"wenty fainter d"warfs (do"wn to V<9) 
"was added in October 1999, follo"wing suggestions that 
metal-enriched stars seemed to be preferentially revealing 
planets (see eg. Laughlin (2000) and references therein). 
These stars all had uvby photometry suggesting metal- 
enrichment. All stars in this fainter sub-sample "were ob- 
served "with a maximum exposure time of 300s regardless 
of observing conditions. As a result, the velocity precisions 
achieved for these targets are not as high as those demon- 
strated for the AAPS main sample (Jones et al. 2002a; 
Butler ct al. 2001). 

2. CHARACTERISTICS OF THE HOST STARS 

HD 30177 (HIP 21850, SAO 233633) is a chromospher- 
ically inactive (logR^^=-5.08; Tinney et al. (2002b)) 
G8V star (Houck & Cowley 1975). HD 30177 was 
observed 105 times by the HIPPARCOS satellite and 
found to be photometrically stable with a standard de- 
viation of 0.017 magnitudes. Its HIPPARCOS par- 
allax of 18.3±0.8mas implies absolute magnitudes of 
Mv=4.72±0.09 (ESA 1997) and Mboi=4.36±0.10 (Cox 
2000) . Stromgren uvby photometry from the General 
Catalogue of Photometric Data-'^^ (Hauck & Mermilliod 
1997), together with the uvby calibrations of (Schuster & 
Nissen 1989) indicate a metallicity of [Fe/H]=+0.20±0.16. 
A detailed metallicity analysis has been performed for all 
the G-dwarfs in the AAPS (Bond 2002), which indicates a 
metallicity of [Fe/H]=-|-0.19±0.09, in agreement with the 
uvby-hased estimate. The b — y colour of HD 30177 indi- 
cates Te//=5320±20K (Hauck & MermiUiod 1997; Olsen 
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1. INTRODUCTION 

The Anglo-Australian Planet Search (AAPS) is a long- 
term planet detection program which aims to perform 
extra-solar planet detection and measurement at the high- 
est possible precision. Together with programmes using 
similar techniques on the Lick 3 m and Keck I 10 m tele- 
scopes (Fischer et al. 2001; Vogt et al. 2000), it provides 
all-sky planet search coverage for inactive F,G,K and M 
dwarfs down to a magnitude limit of V=7.5. Initial results 
from this programme (Tinney et al. 2001; Butler et al. 
2001; Tinney et al. 2002a; Jones et al. 2002a; Butler et 
al. 2002a; Jones et al. 2002b) demonstrate that it achieves 
long-term, systematic velocity precisions of 3ms^^ or bet- 
ter, for suitably stable stars, down to our main sample 
magnitude limit of V<7.5. 

AAPS is being carried out on the 3.9 m Anglo- 
Australian Telescope (AAT), using the University College 
London Echelle Spectrograph (UCLES) and an I2 absorp- 
tion cell. UCLES is operated in its 31 lines mm~^ mode. 
Prior to 2001 September, it was used with a MIT/LL 
2048x4096 15/im pixel CCD, and since then has been used 
with an EEV 2048x4096 13.5^m pixel CCD. Our target 
sample includes 178 FGK stars with 6 < -20° and V<7.5. 
Where age/activity information is available from R^j^ in- 
dices - see for example Henry et al. (1996); Tinney et al. 
(2002b) - we require target stars to have logR'jj^> -4.5 
(corresponding to ages greater than 3Gyr). The observ- 
ing and data processing procedures follow that described 
in Butler et al. (1996) and Butler et al. (2001). 

In addition to our primary sample of V<7.5 dwarfs, a 
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1984). From HD30177's logR^^=-5.08 we estimate in- 
trinsic stellar velocity jitter to be in the range 2-4 ms~^ 
(Saar et al. 1998). 

Figure 1 shows evolutionary tracks at [Fc/H]=-|-0.2 and 
+0.1 for stars of near-Solar mass from the compilation 
of Girardi et al. (2000). The [Fc/H]=+0.1 tracks were 
obtained by interpolation between the [Fc/H]=-|-0.2 and 
[Fe/H]=0.0 tracks in Girardi et al. (2000). Based on these 
models and the metallicity measurements for HD 30177, 
its mass is estimated to be 0.95±0.05Mq, and it would 
appear (in common with HD 73526 and HD 76700) to be 
beginning its evolution off the main sequence. 

HD 73526 (HIP 42282, SAO 220191) is a G6V dwarf 
(Houck 1978) for which no K'f^j^ estimate is currently 
available. It was observed 137 times by HIPPARCOS, 
but its relative faintness means that only a 0.02 magni- 
tude standard deviation limit is placed on its photomet- 
ric stabiHty. Its HIPPARCOS parallax is 10.6±1.0mas, 
which indicates My=4.1±0.2 and Mboi=3.7±0.2 (ESA 
1997; Cox 2000). As for HD 30177 (see above), 
Stromgren uvby photometry indicates a metallicity of 
[Fe/H]=-F0.10±0.16 (Hauck fc Mermilliod 1997; Schus- 
ter & Nissen 1989), while a detailed metallicity analysis 
gives [Fe/H]=+0.11±0.10 (Bond 2002). The GCPD b-y 
colour of HD 73526 indicates To//=5450±20K (Hauck & 
Mermilliod 1997; Olsen 1984). Based on the evolutionary 
tracks shown in Figure 1, HD 73526 is estimated to have a 
mass of 1.05±0.05Mr5. 

HD 76700 (HIP 43686, SAO 250370, LTT3291) is cat- 
alogued as a G6V dwarf by SIMBAD, as G5 by the Henry- 
Draper catalogue (Cannon & Pickering 1918-24) and 
as G8V by HIPPARCOS (ESA 1997). The HIPPAR- 
COS B-V colour (+0.745) would indicate a G8 spectral 
type is more likely than G5 or G6, though the GCPD 
b — y colour (Hauck & Mermilliod 1997) would indi- 
cate a G6 type is more appropriate. No R'^^ esti- 
mate is available. It was observed 109 times by HIP- 
PARCOS, which found it to have constant magnitude 
with a standard deviation of 0.012 magnitude. Its 
HIPPARCOS parallax is 16.8±0.7mas, which indicates 
My=4.3±0.1 andMM=3.9±0.2 (ESA 1997; Lang 1992). 
Stromgren uvby photometry indicates a metallicity of 
[Fe/H]=+0.14±0.16 (Hauck & Mermilliod 1997; Schus- 
ter & Nissen 1989), while a detailed metallicity analy- 
sis gives [Fe/H]=+0.10±0.11 (Bond 2002). The GCPD 
b-y colour indicates a To//=5423±20K (Hauck & Mer- 
milliod 1997; Olsen 1984). Based on the evolution- 
ary tracks shown in Figure 1, HD 76700 is estimated to 
have a mass of l.OO±O.O5M0. There is no published 
log R^jf measurement for HD 76700. Because the veloc- 
ity amplitude observed in this star is small (26±2ms^"'^), 
a Ca HK spectrum was acquired to determine the level at 
which our results could be affected by velocity jitter. Fig- 
ure 2 shows this spectrum for HD 76700 compared to those 
for several objects of similar spectral type, along with their 
spectral types and log R^^ values (Tinney et al. 2002b). 
HD 76700 shows no evidence for a line reversal. Based on 
this comparison we assign an upper limit to the log R^^for 
HD 76700 of -4.9, from which we estimate its intrinsic stel- 
lar velocity jitter to be in the range 3-6 ms^^ or less (Saar 
et al. 1998). 

HD 2039 (HIP 1931, SAO 23205) is a chromospherically 



inactive (R'^^=-4.91; Tinney et al. (2002b)) star clas- 
sified as G2/G3 IV/V by Houck & Cowley (1975). It 
was observed 147 times by HIPPARCOS, but its rela- 
tive faintness means that only a 0.022 magnitude stan- 
dard deviation limit is placed on its photometric stabil- 
ity Its HIPPARCOS parallax is 11.1+1.1 mas, which in- 
dicates My =4.25+0.24 and Mb,,; =4.22+0.25 (ESA 1997; 
Cox 2000). There are multiple measures of this star's 
Stromgren uvby photometry, the mean of which indicate 
a mctalhcity of [Fe/H]=+0.10+0.16 (Hauck & MermiUiod 
1997; Schuster & Nissen 1989). A detailed metallicity 
analysis gives [Fe/H]=+0.10+0.11 (Bond 2002). Theb-y 
colour of HD2039 indicates Te//=5675+20K (Hauck & 
Mermilliod 1997; Olsen 1984). Based on the evolution- 
ary tracks shown in Figure 1, HD 2039 is estimated to have 
a mass of 0.98+0.05Mq. 

3. RADIAL VELOCITY OBSERVATIONS AND ORBITAL 
SOLUTIONS 

Fifteen observations of HD 30177 are listed in Table 1. 

The column labelled "Uncertainty" is the velocity uncer- 
tainty produced by our least-squares fitting procedure. 
This fit simultaneously determines the Doppler shift and 
the spectrograph point-spread function (PSF) for each ob- 
servation made though the iodine cell, given an iodine ab- 
sorption spectrum and an "iodine free" template spectrum 
of the object (Butler et al. 1996). The uncertainty is de- 
rived for each measurement by taking the mean of four 
hundred useful spectral regions (each 2 A long) from each 
exposure. This uncertainty includes the effects of photon- 
counting uncertainties, residual errors in the spectrograph 
PSF model, and variation in the underlying spectrum be- 
tween the template and "iodine" epochs. All velocities are 
measured relative to the zero-point defined by the template 
observation. Only observations where the uncertainty is 
less than twice the median uncertainty are listed. These 
data are shown in Figure 3. The figure shows the best-fit 
Keplerian model for the data, with the resultant orbital 
parameters listed in Table 2. Due to HD30177's faintness 
compared to the AAPS main sample, the residuals about 
the fit (14m s~^) are significantly higher than the baseline 
3m s~^ precision level demonstrated for the main sample 
(Butler et al. 2001; Jones et al. 2002a). 

Eighteen observations of HD 73526 are listed in Table 
3, and they are shown in Figure 4 along with a Keplerian 
fit to the data with the orbital parameters listed in Ta- 
ble 2. Note again the larger residuals (17ms~-'^) compared 
to the AAPS main sample. Twenty-four observations of 
HD 76700 are listed in Table 4, and they are shown in Fig- 
ure 5 along with a Keplerian fit to the data with the orbital 
parameters listed in Table 2. The thirty-six observations 
of HD 2039 are listed in Table 5, and they are shown in 
Figure 5 along with a Keplerian fit to the data with the 
orbital parameters listed in Table 2. 

4. DISCUSSION 

The resultant minimum companion mass for HD 76700 
is Msini = 0.197+0. 017Mjup, with an orbital semi-major 
axis a = 0.049+0. 004 au and eccentricity 6=0.00+0.04. 
This zero eccentricity is consistent with the expectation 
that a planet with a period of just 3.971+0.001 d will 
almost certainly lie in an orbit which has been tidally 
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circularised (Marcy & Butler 1998). The resulting 
orbital parameters for HD 76700 place it amongst the 
planetary companions with the lowest known minimum 
masses. HD 76700, joins HD 49674, HD 16141, HD 168746 
& HD 46375 (Butler ct al. 2G02b; Marcy, Butler & Vogt 
2000; Pepe et al. 2002) in the group of extra-solar plane- 
tary companions with measured minimum masses less than 
a Saturn mass (0.299 Mjup)- 

The remaining three extra-solar planets all have ellipti- 
cal orbits, though the ellipticity of the orbit for the com- 
panion to HD 30177 (e=0.22±0.17) is not different from 
zero with great significance. As further data are acquired 
over the coming years this parameter will become far bet- 
ter constrained, and it is possible that this extra-solar 
planet could turn out to be in a substantially circular or- 
bit. If so this system would join with the other known 
nearly circular systems with gas-giant planets lying in or- 
bits between where the Earth and Jupiter lie in our own 
Solar System (eRet, HD4208, the outer components of 
47UMa, HD 28185 (Jones et al. 2002a)). These extra- 
solar planets would seem to indicate that gas-giants exist 
in nearly circular orbits with semi-major axes all the way 
out to, and beyond that of Jupiter, as confirmed recently 
by the detection of the outer planet in the 55 Cnc system 
(Marcy et al. 2002a). 

Of the twenty metal-enriched stars included as a sub- 
sample along with our main sample in late 1999, five have 
revealed the presence of planetary companions (the four 
planets discussed here - HD 30177, 73526, 76700 and 2039 
- along with the previously known companion to HD 83443 
(Butler et al. 2002a; Mayor et al. 2002)). This gives us a 
lower limit (there may be longer period or lower mass plan- 
ets present which we cannot yet detect) to the discovery 
rate of 25±11% for this "metallicity-biassed" sub-sample. 
This compares with the overall discovery rates estimated 
for the Keck, Lick and AAPS of ~8% (ie 8% of stars sur- 



veyed have planets in orbits within 3.5 a.u. of their host 
stars (Butler et al. 2002c)) - a difference which, while not 
of great statistical significance, is not unexpected given 
that extra-solar planets seem to be being found prefer- 
entially around metal-enriched stars (eg. Reid (2002); 
Laughlin (2000) and references therein). It is also inter- 
esting to note (even if perhaps not statistically significant) 
that all four of the stars in this paper would seem to be 
beginning their evolution off the main sequence (see Fig 

1). 

These results would suggest that the biassing of planet 
surveys toward metal-enriched host stars may offer a ben- 
efit in the planet detection rate. However, such an in- 
creased discovery rate must be balanced against the fact 
that it will produce an inherently biassed sample of extra- 
solar planets. With the total number of extra-solar planets 
still numbering less than a hundred, and the parameters of 
this ensemble of planets still poorly placed in a scheme of 
extra-solar planetary formation and evolution, now is not 
the time for planet searches to begin biassing their large 
surveys in the chase for better "hit rates" at the expense 
of scientific utility. 
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Logio(Teif) Log^(,(T^jj) 

Fig. 1. Evolutionary tracks from the compilation of Girardi ct al. (2000). The Zero Age Main Sequence (ZAMS) is shown as a solid 
line, with evolutionary models at the indicated masses shown as dot-dashed lines. The open circles and uncertainties indicate the stars under 
discussion in this paper, with their estimated uncertainties indicated in brackets. The [Fe/H]=+0.2 tracks axe directly from Girardi et al. 
(2000), while the [Fe/H]=+0.1 tracks were obtained by interpolation between the catalogued [Fe/H]=+0.2 and 0.0 models. 
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Fig. 2. — UCLES spectrum in the region of the Ca H hne for HD 76700 and several comparison objects with log measurements from 
Tinney et ah (2002b). For each object the spectral type and measured logR^^value is shown. Based on these we assign an upper limit to 
the log R'^^value for HD 76700 of -4.9. 
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Fig. 3. — AAT Doppler velocities for HD 30177 from 1998 November to 2002 May. The solid line is a best fit Keplerian with the parameters 
shown in Table 2. The rms of the velocities about the fit is 14ms~^. Assuming 0.95 Mq for the primary, the minimum (Msini) mass of the 
companion is 7.7±1.5Mjup, and the semi-major axis is 2.6±0.9a.u. 
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Fig. 4. — AAT Doppler velocities for HD 73526 from 1999 January to 2002 June. The solid line is a best fit Keplerian with the parameters 
shown in Table 2. The rms of the velocities about the fit is 18 m s~^ . Assuming 1.05 Mq for the primary, the minimum (M sin i) mass of the 
companion is 3.0±0.3Mjup, and the semi-major axis is 0.66±0.05 a.u. 
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Fig. 5. — AAT Doppler velocities for HD 76700 from 1999 February to 2002 June, phased at the best fit Keplerian period of 3.971 d (left 
panel) and plotted unphascd from the last three observing runs in 2002 (right panel). The plotted Keplerian has the parameters shown in 
Table 2. The rms of the velocities about the fit is 6.2ms~^. Assuming I.OMq for the primary, the minimum (Msini) mass of the companion 
is 0.197±0.017Mjup, and the semi-major axis is 0.049±0.004a.u. 
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Fic;. 6. AAT Doppler velocities for HD2039 from 1998 November to 2002 June, with the best fit Keplerian with parameters shown in 
Table 2. The rms of the velocities about the fit is 15 ms~^. Assuming 0.98 Mq for the primary, the minimum (M sini) mass of the companion 
is 5.1±1.7Mjup, and the semi-major axis is 2.2±0.2a.u. 
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Table 1 



Velocities for HD 30177 





RV" 


Uncertainty 


(-2451000) 


(ms-i) 


(ms-i) 


118.0974 


77.9 


10.9 


119.1924 


37.2 


17.1 


121.1514 


49.7 


14.8 


157.1022 


95.4 


8.6 


211.9834 


92.4 


11.3 


212.9660 


92.5 


10.4 


213.9995 


92.3 


9.1 


214.9506 


104.0 


9.3 


525.9973 


26.9 


6.9 


630.9156 


-0.5 


9.9 


768.3296 


-61.5 


10.2 


921.1075 


-117.8 


10.2 


1127.3205 


-168.5 


12.4 


1188.2532 


-172.9 


7.8 


1358.9181 


-173.4 


7.8 



'^Julian Dates (JD) are barycentric. 
Radial Velocities (RV) are barycentric, 
but have an arbitrary zero-point de- 
termined by the radial velocity of the 
template, as described in Section 3 
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Table 2 
Orbital Parameters 



Parameter 


HD 30177 


HD 73526 


HD 76700 


HD 2039 


Orbital period P (d) 


1620±800 


190.5±3.0 


3.971±0.001 


1190±150 


Velocity amp. K (ms~^) 


140±10 


108±8 


25±2 


136±30 


Eccentricity e 


0.22±0.17 


0.34±0.08 


0.00±0.04 


0.69±0.15 




288±40 


207±30 




333±20 


ai sin i (km) 


(3.0±1.5)xl06 


(0.265±0.01)xl0'^ 


1381±1.2 


(1.61±0.45)xl0'^ 


Periastron Time (JD-245000) 


1027±200 


951±12 


1212.9±0.1 


836±150 


Msini (Mjup) 


7.7±1.5 


3.0±0.3 


0.197±0.017 


5.1±1.7 


a (AU) 


2.6±0.9 


0.66±0.05 


0.049±0.004 


2.2±0.2 


RMS about fit (ms"^) 


14 


18 


6.2 


15 



Table 3 
Velocities for HD 73526 





RV° 


Uncertainty 


(-2451000) 


(ms~-^) 


(ms-i) 


212.1302 


23.5 


12.9 


213.1314 


66.8 


13.7 


214.2389 


61.6 


14.9 


236.1465 


48.4 


10.9 


630.0280 


35.5 


11.0 


717.9000 


-120.7 


12.8 


920.1419 


-32.0 


13.4 


984.0378 


78.0 


9.3 


1009.0976 


71.8 


9.7 


1060.8844 


-58.9 


7.4 


1091.8465 


-149.8 


13.6 


1386.9003 


51.8 


6.8 


1387.8921 


64.0 


6.2 


1420.9248 


-14.4 


7.8 


1421.9199 


-9.6 


6.5 


1422.8602 


-9.8 


7.3 


1424.9237 


-19.1 


10.4 


1454.8529 


-109.0 


6.3 



•^As for Table 1 



Four new planets 
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Table 4 

\'el()(ttjes fou lid 7()70() 





RV" 


Uncertainty 


(-2451000) 


(ms-i) 


(ms-i) 


212.1565 


9.2 


6.7 




zLfl 'i 

4:U.O 


7 fi 
1 .u 


214.2583 


-21.3 


8.2 


274.0177 


-18.6 


7.9 


530.1791 


26.7 


8.8 


683.8938 


-19.8 


5.7 


920.1606 


29.7 


7.6 


984.0068 


22.8 


7.1 


1009.0638 


-26.8 


7.5 


1060.9036 


-24.3 


4.7 


1091.8517 


-16.8 


7.4 


1129.8425 


18.2 


6.7 


1359.0760 


-14.3 


6.6 


1386.9145 


-20.0 


3.8 


1387.9170 


23.2 


3.9 


1388.9898 


16.5 


3.4 


1389.8741 


-11.8 


3.9 


1420.9418 


7.3 


4.5 


1421.9331 


-19.3 


3.1 


1422.8793 


-11.0 


2.7 


1424.9439 


0.6 


5.3 


1452.8943 


-4.2 


3.8 


1454.8682 


-3.2 


3.9 


1455.8565 


22.8 


3.1 



As for Table 1 
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Table 5 
Velocities for HD 2039 





RV" 


Uncertainty 


(-2451000) 


(ins~^) 


(ms~^) 


118.0578 


24.0 


8.4 


118.9610 


-23.8 


16.9 


119.9445 


0.0 


11.5 


121.0385 


-33.0 


16.1 


211.9514 


-10.1 


18.7 


212.9234 


3.9 


13.3 


213.9749 


-3.9 


17.0 


214.9171 


-14.1 


13.2 


386.3227 


-54.7 


17.1 


387.2981 


-33.2 


15.3 


411.2293 


-14.2 


15.6 


414.2585 


-16.5 


10.6 


473.0883 


-42.5 


10.6 


525.9286 


-78.7 


13.9 


527.9226 


-28.9 


11.0 


745.2702 


-51.8 


17.4 








828.9944 


-49.3 


10.4 


829.9757 


-45.1 


11.7 


856.0702 


-33.5 


14.6 


919.9434 


-25.7 


14.7 


920.9672 


-22.5 


14.5 


1093.2947 


136.2 


UA 


1127.2341 


98.9 


15.1 


1151.2230 


65.7 


8.4 


1152.0860 


66.4 


9.2 


1154.2124 


74.6 


11.4 


1187.0957 


55.1 


9.7 


1188.0300 


48.7 


8.4 


1189.1502 


58.6 


11.6 


1190.0932 


48.5 


8.4 


1422.3281 


-9.7 


11.2 


1425.3322 


-3.2 


7.7 


1455.2853 


-29.4 


6.1 



As for Table 1 



